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Abstract 

For given matrix A £ Z dxn , the set Pa,6 — { z '■ Az = b,z £ describes the preimage or 
fiber of b £ Z d under the Z-linear map : Z™ — > 7h d , a; Ax. The fiber P\ b is called atomic, 
if Pk,b — ^A.bi + ^A.b 2 implies b = bi or 6 = 62- In this paper we present a novel algorithm to 
compute such atomic fibers. An algorithmic solution to appearing subproblems, application 
to integer programming, and computational examples are included as well. 

1 Introduction 

Decomposition of rational polyhedra is at the heart of several interesting applications. However, 
there are different definitions of decomposability depending on the application. These definitions 
mainly differ in the treatment of the (integer) points of the polyhedron. 

The simplest notion is that of linear decomposition of polyhedra. Two polyhedra P,Q C R™ are 
called homothetic if P = XQ + t for some A > and t £ W l . Here, a polyhedron P is called 
indecomposable, if any decomposition P = Qi + Q2 implies that both Q± and Q2 are homothetic 
to P. It can be shown that there are only finitely many indecomposable rational polyhedra that 
are not homothetic to each other. For further details on this type of decomposition we refer the 
reader for example to □ US EH H3 

Let us now come to a bit more restrictive decomposition. Here we consider only polyhedra of the 
form {x € W 1 : Ax < b} for given matrix A G Z dxn and varying b € Z d . To emphasize that we only 
consider integer right-hand sides, we say that a polyhedron P is integrally indecomposable, if any 
decomposition P — Qi + Q2 (into polyhedra with integer right-hand sides) implies that both Q\ 
and Q2 are homothetic to P. This decomposition is more restrictive than the linear decomposition, 
since only such polyhedra Qi and Q2 are allowed that have an integer right-hand side. In 0, 
the authors showed that properties like TDI-ness of each member of a family of systems Ax < b, 
b € Z d , can be concluded from TDI-ness of the integrally indecomposable systems. 

Another important application of integral decomposition of polyhedra is that of factorizing a 
multivariate polynomial. Here, one considers only polyhedra of the form {x G K" : Ax < b} for 
given matrix A £ °z dxn and varying b £ 7L d , where each polyhedron is integer, that is, where each 



polyhedron has only integer vertices. The reason for this restriction is the simple observation that 
the so-called Newton polytope Newt(/) := conv(a : x a G supp(/)) associated to a polynomial / is 
integer by definition. Moreover, the relation / = gh among three polynomials /, g, and h implies 
Newt(/) = Newt(g) + Newt(/i), a theorem due to Ostrowski. For more details on this subject see 
for example |12| and the references therein. 

In this paper, we will consider another notion of decomposability. In contrast to considering all 
points of a polyhedron, we restrict out attention to its set of lattice points. More formally, we 
consider family of sets of lattice points 

P Ab := {z : Az = b,z G Z+}, 

where the matrix A G Z dxn is kept fix and the right-hand side vector b G Z d varies. Note that 
P{ b is exactly the preimage or fiber of b under the linear map Ja ■ Z" — > Z d , x i— > Ax. We call a 
fiber P A b indecomposable or atomic, if PjJ 6 = P l A bi + P A (, 2 implies b — b\ or 6 = &2 ■ Note that 
P4 6 = P^ bi + P A j, 2 means that every lattice point of P£ & is the sum of a lattice point of P£ bi 
and a lattice point of P4 62 . This is indeed a very strong condition, but again it can be shown that 
there are only finitely many (nonempty) atomic fibers for a given matrix A 9 . Atomic fibers were 
used in pP to construct strong SAGBI-bases and receive more recent attention in the computation 
of minimal vanishing sums of roots of unity |14|. Note that atomic fibers are not only minimal 
(with respect to decomposability) within the given family, but also generate every fiber P l A b in 
this family as a Minkowski sum P A b = ©^ =1 otiP A b . , 014 E Z+. 

Below, we present a new algorithm to compute atomic fibers. In fact, this algorithm can also 
compute - what we call - extended atomic fibers. For this, we call 

Q A>b :={z:Az = b,ze Z"} 

an extended (integer) fiber and we call it atomic, if (Q A b n Oj) = (Q A b tl Oj) + (Q A b2 D Oj) 
holds for all the 2" orthants Oj of K™, then b = 61 or b = 62- Also here it can be shown that there 
are only finitely many (nonempty) extended atomic fibers for a given matrix. Also this very strong 
notion of decomposability has an application: the set Hoo constructed in [H] for use in two-stage 
stochastic integer programming is in fact the set of extended atomic fibers of the family of extended 
fibers 

{x, y : x = b, Tx + Wy = 0, x G Z m , y G Z"} 
where T and W are kept fix and where b varies. 

One algorithm to compute atomic fibers via certain standard pairs was presented in £Q. Via this 
algorithm, the authors computed the atomic fibers of the twisted cubic by hand(!), see Example 
12.31 Recently, the computation of atomic fibers appears a sub-problem or are being investigated in 
the design of telecommunication networks with uncertain demand and capacity. 

Before we present our algorithm, we include in Section|3|an application of atomic fibers: the integer 
program min{c T 2: : z G P A b } can be solved via the optimal solutions of the integer programs whose 
right-hand sides are those of the atomic fibers. 



2 



2 Atomic Fibers 



Let us now start our treatment with a formal definition of atomic fibers. 

Definition 2.1 For u, v G R n we say that u C v ifu^'v^) > and \u^ \ < \v^\ for all compo- 
nents j = 1, . . . , n, that is, u belongs to the same orthant as v and its components are not greater 
in absolute value than the corresponding components of v. 

For S,T,U C R n we say that S = T © U if for every z e S there are z\ G T and Z2 £ U with 
Zi,%2 and z — z\ + Z2- 

For A G Z dXTl and b e Z d we call P{ b := {z : Az = b, z G Z^} a fiber and Q A b := {z : Az = 
b, z G Z"} an extended fiber o/A. 

VFe call P A b atomic, i/P| b cannot be written as P A b — P A h ®P A b2 with two other fibers P A b 
and P A,b 2 ■ 

Analogously, we call Q A b atomic, if Q A b cannot be written as Q A b = Q A bi © Q A &2 with two 
other extended fibers Q A b and Q A b<2 . 

By F{A) and by E{A) denote the atomic and the extended atomic fibers of A, respectively. 

Note that for every integer matrix A there are at most finitely many (extended) atomic fibers. This 
fact is based on the following nice theorem. 

Theorem 2.2 (Maclagan, |9j) Let T be an infinite family of monomial ideals in a polynomial 
ring k[x\, . . . , x n \. Then there must exist ideals I, J G X with I C J . 

To apply this theorem in our situation of atomic fibers, associate to the fiber Q A b the monomial 
ideal 2^,6 ■= (x^ z ' z ' : Az = b, z G Z") C Q[xi, . . . , X2 n ]. Q A b is an extended atomic fiber if and 
only if 7a, b is inclusion-maximal among all these ideals. 

For the treatment below we like to point out here that the relations P A b — Q A b (~l R™ and 
P\ bi ® Pa 62 = P A bi + Pa b 2 — Pa b x +b 2 are eas y to ver ify- Another simple fact to notice is that 
QAM+b 2 = Qam + Q'am whenever Q Abl # or Q AM £ 0. 

Example 2.3 In ^j, it was shown how atomic fibers could be used to construct strong SAGBI 
bases for monomial subalgebra over principal ideal domains. As an example, they computed the 
atomic fibers of the matrix A = (0123) hand via an approach different from the one we present 
below. 

In the table below, we list the right-hand sides and all (finitely many) elements in these 18 atomic 
fibers. (Using 4ti2 OV ; in which the algorithm presented below that computes (extended) atomic 
fibers has been implemented, we can not only verify this result, but also quickly find the 51 extended 
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fibers associated to this matrix.) 
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(1,0, 3,1), (0,2, 2,1), (0,1, 4,0)} 
(2, 0,3,0), (1,2, 2,0), (0,4, 1,0)} 
(1,0, 3, 2), (0,2, 2, 2), (0,1, 4,1), (0,0, 6,0)} 
(3, 0,3,0), (2, 2, 2,0), (0,6, 0,0), (1,4, 1,0)} 



Thus, for example, the fiber given by the right-hand side 
decomposed into atomic fibers as 



, 7) is not atomic, since it can be 
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This can be quickly verified by looking at the elements in these fibers: 

{(2, 1,0,2), (2,0,2,1), (1,1,3,0), (1,2, 1,1), (0,4, 0,1), (0,3, 2,0)} 
= {(0,1,0,1), (0,0,2,0)} ©{(2,0, 0,1), (1,1, 1,0), (0,3, 0,0)}. 

Indeed, we have 

(2,1,0,2) = (0,1, 0,1) + (2, 0,0,1), 
(2,0,2,1) = (0,0, 2,0) + (2, 0,0,1), 

(1.1.3.0) = (0,0, 2,0) + (1,1, 1,0), 

(1.2.1.1) = (0,1, 0,1) + (1,1, 1,0), 
(0,4,0,1) = (0,1, 0,1) + (0,3, 0,0), 
(0,3,2,0) = (0,0, 2,0) + (0,3, 0,0). 

In Example l2 . 31 above . it was easy to verify whether a given fiber is a summand in the decomposition 
of another fiber by simply checking the finitely many elements in the fiber for a decomposition. If 
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the fibers are not bounded, however, this would not give a finite procedure. The following Lemma 
tells us how to solve this problem via the (finitely many!) C-minimal elements in the given fibers. 

Lemma 2.4 Let Q AM ^ and Q Ab2 ^ 0- Then Q Abi+b2 = Q AMl © Q A .b 2 and onl y f or ever y 
^-minimal vector v G Q I A.b 1 +b 2 there is a vector w G Q A bl with ti)C«, 

Proof. Let v G Q I Ab 1 +b 2 ari< ^ ^ ^ e Q I A.b 1 +b 2 ^ e ^-minimal in Q Abl+b2 with v C v. Thus, 
by the assumption in the lemma, there is some w G Q A bi such that w \Z v. Clearly, this yields 
v — w G &2 and w — w C v. 

We now claim that v = (w + v — v) + (v — w) with w + v — v G Q A bi , v — w E Q A h2 , w + v — v Qv, 
and ti - w C u, is a desired representation of v. The first two relations are trivial, if we keep in 
mind that Av = Av = b, Aw = b\, and b = b\ + b 2 - We get the other two relations as follows: 

(a) w + v — v ^ v + v — v — v, since by construction w and v — v lie in the same orthant, and 

(b) v — w C v C v, since w C v. 

Thus, we have constructed for arbitrary v € Q A bl+b2 a valid representation of v as a sum of two 
elements from Q A b and Q A b2 that lie in the same orthant as v. This concludes the proof. □ 

Using this lemma repeatedly, we are now able to find, for a given right-hand side b, a decomposition 
Qa b = a iQ l A b ■■> a i ^ ^+ ' that is, we can find a decomposition of a fiber into a sum of 

atomic fibers. Note that replacing Q by P and E(A) by F(A) in the algorithm below would yield 
a decomposition of the atomic fiber: P A b = (J)^! ctiP A b . , a, G Z+. 

Algorithm 2.5 (Algorithm to decompose extended fiber into extended atomic fibers) 
Input: A, E(A) = {bi, . . . , b k } 

k 

Output: ai,..., afe such that Q A b — @ aiQ^ b . 

' i=i 

oti := . . . := a fe := 
for i from 1 to k do 

while Q A ,b = Qam ® QA,b-bi 4q 
b:=b-h 
ai := a t + 1 
return ot\, . . . , a k - 

It remains to state an algorithm that computes the finitely many C-minimal elements in P A b and 
in Q A 6 , respectively. One simple way to compute these elements is via (truncated) Hilbert and 
Graver bases. For this we note that the C-minimal elements in P A b correspond to those elements 
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in the Hilbert basis of the cone {(z,u) : Az — bu = 0, (z,u) G that have u — 1, while 

the C-minimal elements in Q A b correspond to those elements (z,u) in the Graver basis of the 
matrix (A\ — b) that have u = 1. Finiteness of such Hilbert and Graver bases for rational data 
A, b immediately implies that there are indeed only finitely many C-minimal elements in P A b and 
Q A b , respectively. This, however, could have been already concluded from the fact that C-minimal 
elements in Q A b are in one-to-one correspondence with the (finitely many) minimal generators of 
the monomial ideal lA,b- 

3 Solving Integer Programs via Atomic Fibers 

Algorithm 12.51 puts us in the position to solve integer programs with the help of atomic integer 
programs whose right-hand sides correspond to those of the atomic fibers. Therefore, these atomic 
programs encode already the solutions to a whole family of integer programs (for changing right- 
hand side). 

Definition 3.1 (Atomic Programs) If P l A b is an atomic fiber of A G 1, dxn , then we call 

min{c T z : z G P A b } = min{c T z : Az = b, z G Z™ } 
the atomic (integer) program associated to P Ab - 

Once we know optimal solutions to these finitely many atomic programs, we can solve the integer 
program for any given right-hand side b by simply decomposing the fiber P A b into a sum of atomic 
fibers. 

Lemma 3.2 Given A G Z dxn , c G R n , and F(A) = {bi,...,b k }, let z° pt , . . . , z° k pt be optimal 
solutions to 

min{c T z : Az = bi, z G Z™ }, i = 1, . . . , k, 
and consider any solvable problem 

min{c T z : Az = b, z G Z n }. (1) 
V -^A b = ©i=i a «^A 6 > a i 6 then z :— $^ i=1 ctiZ° pt is an optimal solution to f7|). 
Proof. First we show that z is a feasible solution of (JTJ- Clearly, z — ^ i=1 ctiZ° pt > 0. Moreover, 

Az = A ( aiZ?*) =J2 a i ( Az T) = E a ^ = b - 

\i=l / i=l i=l 

Therefore, z is a feasible solution of 0J. 



G 



Let z be any feasible solution of Q. From the decomposition b = © i=1 Qfi-P^ b . with G Z + 
we conclude that there are Zi G P4 b . , i = 1, . . . , k, such that z = Y]j—i ctiZi- But now we have 



1 / z— 1 i— 1 \i— 1 



Therefore, z is optimal for □ 



4 Computation of Atomic Fibers 

In the following we show how to compute the finitely many (extended) atomic fibers of a matrix. 
We use the algorithmic pattern of a completion procedure. By e\, . . . , e n we denote the unit vectors 
in R n . 

Algorithm 4.1 (Algorithm to compute atomic fibers) 
Input: F := {Aei, . . . , Ae n , —Ae\, . . . , —Ae n } 

Output: a set G, such that {Q J A b : b G G} contains all extended atomic fibers of A 
G:=F 

C := U {f + g} (forming S-vectorsJ 

f,g£G 

while C ^ do 

s := an element in C 

C:=C\{s} 

f := normalForm(s, G) 

iff =£0 then 

G:=GU{/} 

C := C U 1J {/ + g} (adding S- vectors,) 

gee 

G:=GU{0} 
return G. 

Behind the function normalForm(s, G) there is the following algorithm. 

Algorithm 4.2 (Normal form algorithm) 
Input: s, G 

Output: a normal form of s with respect to G 
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while there is some g € G such that Q A 

s := s — g 
return s 



Lemma 4.3 Algorithm \4-l\ terminates and computes a set G such that {Q A t '■ & G G} contains 
all atomic fibers of A. 

Proof. Associate with b the monomial ideal I am '■= {x^- z+ ' z ' : Az — b,z 6 Z n ) C Q[xi, . . . , ar2n]- 
Algorithm 14.11 generates a sequence {fi,f%, ■ ■ ■} in G \ F such that ^. ^ ^ © /■-/< 
whenever i < j. Thus, the corresponding sequence {IaJi , Iaj 2 i •• ■} °f monomial ideals satisfies 
I(A,fj) % I(A,fi) whenever i < j. We conclude, by Maclagan's theorem j^J, Theorem 12.21 above . 
that this sequence of monomial ideals must be finite and thus, Algorithm 14. II must terminate. 

It remains to prove correctness. For this, let G denote the set that is returned by Algorithm 14.11 
Moreover, let Q A r be an extended atomic fiber of A with 6^0. We will show that b £ G. 

Since F C G\ {0}, we know that Q Ab = YlQAb f° r finitely many (not necessarily distinct) 
bj € G \ {0}. This implies in particular, that every z £ Q A b can be written as a sum z = J2 v j 
with Vj £ Q A b _ . We will show that we can find bj £ G such that every z £ Q A b can be written 
as a sum z = J2 v j with Vj £ ^ C Q A ^ and with Vj Q z. This implies 5 = (B ^ ■ Since 
Q A t is atomic, and thus indecomposable, this representation must be trivial, that is, it has to be 
Q A b = Q A 7, and therefore we conclude b £ G. 

For A £ Z dxn and b £ Z rf let 6 denote the set of C-minimal elements in Q A b . It is sufficient 
to find Q A b . £ G such that every z £ b can be written as a sum z = Vi with Uj £ b . , 
where all i?i lie in the same orthant as z (as this then implies Mj Cz), To see this, take any non- 
minimal z' £ Q A t \ M A 7 and let z £ M A b be such that z C. z' . Then, clearly, h := z' — z £ Q A 
and every representation z = m with Uj £ b . and «j C z readily extends to a desired valid 
representation z' = z + h = (v- Lo + h) + Xa^j v ii since Vi + h £ Q A b , and since Ui and h and 
thus also Vi + h lie in the same orthant as z and z' . Consequently, h C z', as required. 

By the Gordan-Dickson Lemma, M A b is finite. Thus, let denote its elements by M A b = {z\, . . . , Zk}- 
From all representations Q A b = Y^jeJ Qa bj w hh bj £ G \ {0} and where Zj = Y^jeJ v i,j an d 
£ Q A b ., i — 1, . . . , k, choose one such that the sum 

EEKjII* ( 2 ) 

is minimal. By the triangle inequality we have that 

k k 

EE ii^-ii^ Ewii- ( 3 ) 

i=l j"e,7 i=l 
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Herein, equality holds if and only if all vtj have the same sign pattern as Zi, i — 1, . . . , k, that is, if 
and only if we have Vij C Zi for all i and all j. Thus, if we have equality in J2J| for such a minimal 
representation Q A b = Qa b' then v i,j G fe. an d ^ij E Zi for all occurring v^j, and we 

are done. 

(It should be noted that we have required bj G G \ {0} for all appearing bj, that is in particular, 
bj 7^ 0. Those bj will be sufficient to generate all [--minimal elements in the extended fiber Q A b . 
The remaining elements in Q A g we get by adding elements from Q A as mentioned above.) 

Therefore, let us assume that 

k k 

E£lkilli>X> w 

i=l j£.J i=l 

In the following we construct a new representation Q Ab = Y^j'ej' Qav wnose corresponding 
sum (J2J) is smaller than the minimally chosen sum. This contradiction proves that we have indeed 
equality in © and our claim is proved. 

From J3J we conclude that there are indices ioiiijj2 and a component m S {1, . . . , n} such that 
v iTji ' v i™h < 0- The extended fiber Q A ^ +b j2 was reduced to Q A by sets Q A b n , j" G J". This 
gives representations 

v i.ji + u i,i2 = ^ ' w i,j" i ^i.i" G Q^.b ■// ' E v i,ji + u iii2 

for i = 1, . . . , k. As all he in the same orthant of R" as Ujjj + Vij 2 , we get 

II £ ^J" II 1 = Ikiji +^J2 111 < 11^,^111+ ||Wi,i 2 ||l, 

j"ej" 

with strict inequality for i = iq. 

Thus, replacing in Q Ab = £\ £J Qa^ the term QA,b n + Qa,^, h y Y,j» e j» Qa,^,, > we arrive at a 
new representation g = Ylj'eJ' Qa b-t whose corresponding sum J2J) is at most 

k k 

EE \Krh< EE ii^-iin 

contradicting the minimality of the representation Q A g = X^jeJ b ■ This concludes the proof. 

□ 

Having an algorithm available that computes all extended atomic fibers, we can of course use it to 
compute atomic fibers: If P l A b is atomic then so is Q A b , as any decomposition of Q A b , restricted to 
the nonnegative orthant, would give a decomposition of P{ b . This way of computing atomic fibers, 
however, is pretty slow, since there are far more extended atomic fibers as atomic fibers. A similar 
behavior can be observed when one extracts the Hilbert basis of the cone {x : Ax = 0, x G K™} 
from the Graver basis of A, as the Graver basis is usually much bigger than the Hilbert basis one 
is interested in. 
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In jS] it was shown that one can reduce this difference in sizes by considering only a subset of k 
variables, by computing all indecomposable vectors in the set {x : Ax — 0, x 6 Z+ x Z n_fe }, and 
by iteratively adding a new variable. The main observation exploited is that the indecomposable 
vectors in {x : Ax = 0,iG Z+ x Z n_fe } are nonnegative on the first k components and in addition 
generate all indecomposable vectors in {x : Ax = 0,x e Z^ +1 x Z Tl ~ fe_1 } as a nonnegative integer 
linear combination. This keeps the sets from which the indecomposable vectors are extracted 
comparatively small and leads to a tremendous speed-up in the algorithm. With this algorithm, 
bigger Hilbert bases, even with more than 500,000 elements, can be computed nowadays. 

More interesting for our problem of finding atomic fibers is the fact that a similar "project-and- 
lift" idea can be applied in our situation here. Consider the set of "partially extended" fibers 
Qa b k := { z : ^ z — b,z <E 1>\ x Z"~ fe }, where A is kept fix and where b varies, and let us 
restrict the conditions for decomposability to the first k components. Now we use Algorithm 14.11 
to compute a set of indecomposable/minimal partially extended fibers Q A b± k , . . . , Q A b k . They 
have the property of being nonnegative on the first k components in each vector and of generating 
every other partially extended fiber Q A b k as a (nonnegative integer) Minkowski sum. This latter 
fact can be used in the computation of indecomposable extended fibers Q A b k+1 by defining the 
input set to Algorithm 14. II to be F := {b\, . . . , b s , Aef--i, ■ ■ ■ , Ae n , ~Aek+i, ■ ■ ■ , —Ae n }. Note that 
we do not need to add — b\, . . . , — b s here due to the nice generating property that we have already 
achieved for Q^ jbljfc) . • . , Q A ,b sl k- 

Clearly, when k = n we have Q A bn = P l A b and we have found the atomic fibers of A. Let us now 
demonstrate the capability of this algorithm by considering the following example. For this note 
that both Algorithm 14. II and the "project-and-lift" version of it to compute atomic fibers directly 
are implemented in 4ti2. 

Example 4.4 One example that appears and was solved in \14\ is the computation of the atomic 
fibers of the matrix 

( 1 -1 -1 1 \ 
1-1 0-1 1 
1 -1 -1 1 
V 1 -1 -1 1/ 

This matrix corresponds to a certain problem on 3 x 3 tables and has in fact 31 atomic fibers and 
79 extended atomic fibers, as can be easily verified with 4ti2. 
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The next higher problem on 4 x 4 tables leads to the matrix 
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4ti2 computes 12, 675 atomic fibers for this matrix. 



5 Conclusions 

In this paper we presented a completion procedure to compute atomic and (extended) atomic 
fibers of the form P Ab = {z : Az = b, z £ Z" } and Q A b — {z : Az — b, z £ Z™}, respectively, and 
demonstrated its ability to solve non-trivial examples. 

If the number n of variables is not fixed, the problem of computing atomic fibers is NP-hard, since 
it is at least as hard as computing Hilbert bases of rational cones, which are already exponential in 
size. However, as Barvinok and Woods showed in the Hilbert basis (and thus also the Graver 
basis) can be encoded via a short (polynomial size) rational generating function if the dimension is 
kept fixed. This implies that the elements in each particular fiber P A b = {z : Az — b, z £ Z" } (and 
equally in each extended fiber) can also be computed in polynomial time, if the dimension is kept 
fix. It is an open question whether a similar short rational function expression for the (right-hand 
sides of the) atomic fibers exists and if it is polynomial time computable in fixed dimension, that 
is, we ask for a short rational function that encodes the generating functions 



E zb and E 

b£F(A) beE(A) 



z b . 



At this point we would like to mention that one may of course use Algorithm 14.11 also for the 
problem of finding indecomposable extended fibers (and thus also indecomposable fibers) among 
the elements in the family of extended fibers Q A b — {z : Az — b, z £ Z"} where A is kept fix 
and where b is allowed to vary only on a sublattice A of 1 d . One only needs to replace the input 
set F = {Ae\, . . . , Ae n , —Aei, . . . , —Ae n } by a symmetric generating set F — {db&i, ■ ■ ■ , ±&fe} of A 
with the additional property that Q A b — {z : Az = b, z £ Z n } is not empty for all b £ F . 

Acknowledgment. The author wishes to thank Jesus De Loera for some useful suggestions on 
improving the presentation of this paper. 
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